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Abstract 

We present new sets of fragmentation functions for neutral kaons, both at leading 
and next-to-leading order. They are fitted to data on inclusive K° production in 
e+e- annihilation taken by MARK II at PEP (y/s = 29 GeV) and by ALEPH at 
LEP. Our fragmentation functions lead to a good description of other e + e~ data 
on inclusive K° production at various energies. They also nicely agree with the K s 
transverse-momentum spectra measured by HI at the DESY ep collider HERA, by 
UA5 at the CERN SppS Collider, and by CDF at the Fermilab Tevatron. 
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1 Introduction 



Recently, precise data on inclusive 7r ± , K ± , and unspecified-charged-hadron production 
in e + e~ annihilation at the Z resonance has been published. Using this new data and 
similar data from lower centre-of-mass (CM) energy (-^/s = 29 GeV), we constructed new 
sets of fragmentation functions (FF) for charged pions and kaons at leading order (LO) 
and next-to-leading order (NLO) [|IJ. These new parameterizations were tested against 
data on tt ± , K^, and charged-hadron production in e + e~ annihilation at various energies 
and data on single-charged-hadron production in small-Q 2 ep scattering at HERA, which 
presents a nontrivial check of the factorization theorem of the QCD-improved parton 
model. 

Besides charged pions and kaons or just charged hadrons, Kg mesons are easily de- 
tected through their dominant decay into 7r + 7r~ pairs. The ALEPH [fj, DELPHI ||, 
OPAL |4j] , and L3 || Collaborations at LEP have recently reported on their high-statistics 
analyses of inclusive single K° production.^ 

Following our strategy of constructing FF for charged pions and kaons, we shall com- 
bine this new data on K° production at the Z resonance with the rather precise data 
taken at y/s = 29 GeV by the MARK II Collaboration [§ at PEP to obtain FF for the 
neutral kaons. Owing to the factorization theorem, the same FF can be used to predict 
cross sections of inclusive single K° production at high transverse momenta (pr) in other 
processes like ep and pp scattering. The functions characterizing the fragmentation of 
gluons, u, d, s, c, and b quarks (antiquarks) into neutral kaons contribute quite differently 
in these processes as compared to e + e~ annihilation. For example, in e + e~ annihilation 
at the Z resonance, all five quarks are directly produced, whereas the gluon does not 
directly couple to the electroweak currents. The gluon only contributes in higher orders 
and mixes with the quarks through the Q 2 evolution. On the other hand, in the case of 
inclusive light-meson production at moderate pt in high-energy pp collisions, the cross 
section is dominated by gluon fragmentation jTj . In ep collisions with almost real photons 
at HERA, the situation is mixed. In the lower pt range (pr;^15 GeV), inclusive sin- 
gle hadron production proceeds dominantly via the resolved photoproduction processes 
99 ~^ 99i Q9 ~ * 9Q-I an d qg —* qg, where the first and second partons originate from 
the virtual photon and the proton, respectively, while the third one fragments into the 
outgoing hadron ||. Direct photoproduction only plays a significant role at larger p^ |§. 
Therefore, quark and gluon fragmentation should give comparable contributions even at 
small pt- 

In our previous work on FF for charged pions and kaons M, we could exploit the 
information from tagged three-jet events in e + e~ annihilation to constrain the gluon 
fragmentation into charged hadrons, which constrained also that into charged pions and 
kaons. Unfortunately, such information is not yet available for inclusive K° production in 
e + e _ annihilation. Thus, we shall have to resort to the information on gluon fragmentation 
into charged kaons which we extracted in Ref . Q . 

1 Unless stated otherwise, we shall collectively use the symbol K° for the sum of Kg and if£ (or K° 
and K°). 
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Another problem that requires special attention is related to the distinction of the 
different quark flavours in K° fragmentation. In our recent analysis of and K frag- 
mentation |T[ , we had some information on the fragmentation of specific flavours at our 
disposal. Preliminary measurements of charged-hadron production by the ALEPH Col- 
laboration |IIJ distinguished between three cases, namely the fragmentation of (i) u, d, 
s quarks, (ii) b quarks only, and (iii) all five flavours (u, d, s, c, and b). This enabled us 
to remove the assumption that the s, c, and b quarks fragment into charged pions (kaons) 



in the same way, which we had made in our earlier work Although equivalent in- 
formation is still lacking for K° fragmentation in e + e~ annihilation, we shall follow the 
approach of our recent work on 7r ± and K production [|l[] , where no additional identities 
between the FF of different quark flavours were imposed, except those following from the 
flavour content of the produced mesons. Should it turn out that the relative importance 
of the different flavours cannot yet be pinned down so reliably, then this will not be due 
to a shortcoming of this specific procedure; this would just signal that more detailed data 
is indispensable in order to determine the differences in flavour of the FF more accurately, 
leaving room for further improvements. 

It is the purpose of this work to make use of the new K° data by ALEPH [[| together 
with the K° data by MARK II f| to construct new LO and NLO sets of FF, only 
identifying the FF of the d and s quarks and imposing no constraints on the quarks 
otherwise. At the starting scale, Qo, we shall take the gluon FF of the neutral kaons to be 
equal to their charged counterparts. The recent data from DELPHI 0, OPAL M, and L3 
H agree with the ALEPH data and will not be used in our fit. A comparison of all four 
data sets may be found in a report by OPAL ||. We choose the ALEPH data, since, in 
the region of relatively large x, which we are mainly interested in, it has a slightly smaller 
total error than the data from DELPHI and OPAL. The data from L3 does not extend 
to x values in excess of 0.24 and is thus less useful for our purposes. 

Our new K° FF sets will be tested against older data from e + e~ colliders with lower 
CM energies. Furthermore, we shall calculate the px distributions of K s mesons produced 
inclusively in ep and pp collisions at various CM energies and compare them with prelim- 



inary HI data |T2| and with data from UA5 JTSJ and CDF in order to check whether 
the gluon fragmentation and the relative importance of the various quark flavours are 
realistically described. 

An alternative way of constructing FF is to fit to data generated by well-established 
Monte Carlo event generators rather than experimental data. This avenue has just re- 
cently been taken in Ref. [15 , where, among other things, a NLO set of Kg FF has been 
presented. This offers us yet another opportunity to test our K° FF, namely against 
Monte Carlo output. We shall report the outcome of such a comparison later on. 

The LO and NLO formalisms for extracting FF from e + e~ data are comprehensively 
described in our previous works [|T],[TT|[ and will not be reviewed here. Also, the formulae 
that are needed to calculate the cross sections of inclusive single hadron production in 
ep collisions (with almost real photons) and in pp collisions may be found in earlier 
publications |1]|||J. The NLO formulae in these references are based on the works by 
Aversa et al. ]n| (resolved photoproduction and pp collisions), Aurenche et al. [[T7|] 
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(direct photoproduction), and Altarelli et al. JT8|] (e + e~ collisions). 

This paper is organized as follows. In Sect. 2, we shall describe the actual analysis 
and present our results for the K° FF. We shall also check these FF against e + e~ data 
at lower energies which we did not use in our fits. Furthermore, we shall compare the 
calculated inclusive K s cross sections for ep and pp collisions with HI, UA5, and CDF 
experimental results. Our conclusions will be summarized in Sect. 3. In the Appendix, 
we shall list simple parameterizations of our FF sets for inclusive K° production. 

2 Results 

For our analysis, we select the data on K° production taken at energy y/s = 29 GeV by 
the MARK II Collaboration at PEP and those collected at y/s = M z by the ALEPH 
Collaboration at LEP M. This data comes in the form (l/a^^da/dx as a function 
of x = 2E K o/y/s, where y/s and E K o are the e + e~ and K° energies in the CM system, 
respectively. The data from MARK II and ALEPH lie within the ranges 0.036 < x < 0.69 
and 0.003698 < x < 0.8187, respectively. These and other e + e~ experiments present 
inclusive cross sections for K s + K° L (or K° + K ), i.e., the sum of the two individual 
rates. We adopt this convention, i.e., our FF refer to the fragmentation of any given 
parton into K s and (or K° and K°). For the fitting procedure, we use the x bins 
in the interval between x min = max(0.1, 2 GeV / 'y/s) and x max = 0.8 and integrate the 
theoretical functions over the bin widths, which is equivalent to the experimental binning 
procedure. The restriction at small x is to exclude events in the non-perturbative region, 
where mass effects are important. Very-large-x data suffer from huge uncertainties, so 
we prefer to disregard the few data-points above x max . As usual, we parameterize the 
x dependence of the FF at the starting scale Q as 

Df (x,Q 2 ) = Nx a (l-x)P, (1) 

where a stands for any quark flavour or the gluon. We impose the condition 

Df°+^(x,Q 2 ) = Df + *°(x,Q 2 ). For all the other quark FF, we take N, a, and 
to be independent fit parameters. 

As mentioned above, the e + e~ data on inclusive single particle production does not 
well constrain the gluon FF, which, however, plays an important role in ep reactions and 
even more so in pp processes. Since, at present, there exists no additional information on 
gluon fragmentation to neutral kaons in e + e~ annihilation, we fall back on the results on 
the fragmentation of gluons into charged kaons obtained in our recent analysis [Q]. We 
argue that the supposedly flavour-blind gluon should fragment into charged and neutral 
kaons at the same rate, and identify the corresponding FF. Later on in this section, we 
shall demonstrate in more detail that, in want of better data, this is a sensible prescription. 

Of course, the data on 7r ± and production have much better statistics than the K° 
production data under investigation in this paper. For this reason, and for compatibility 
with our 7r^~ and K ± sets, we do not fit A^jg anew, but adopt the values determined in 
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Ref. [@, A|| = 107 MeV (195 MeV) in LO (NLO). We are thus left with 12 independent 
fit parameters. 

The quality of the fit is measured in terms of the Xd.o.f. f° r & U selected data points. 
The technical procedure to determine these 12 parameters, using well-tested numerical 
techniques of multidimensional optimization [BJ, is similar to our earlier work ||T|. As in 
Ref. 0, we choose Q = \[2 GeV for the u, d, and s quarks, Q = m(r] c ) = 2.9788 GeV 
for the c quark, and Qq = m(T) = 9.46037 GeV [^OJ for the b quark. Our results are 



listed below. For the sum of K° and K° , we find 

£>f 0+W °' hO \x,Ql) = 0.51a;- - 841 (1 - a;) 1 - 55 , 
Df° + x°< LO \x,Ql)=Di K0+lf ' hO \x,Q 2 ) = 1.47x-°- 691 (l 
Df 0+lf ^°\x,Ql) = 1.00 a;- - 738 (1 -a;) 2 ' 93 , 
Df ° +lf ' h0 \x, Ql) = 0.68 x-°- 598 (1 - x) 1 93 , 



D (K0+K .LO)^ q2 } = .43 3.-0.374 (1 _ x) 2.69 (2) 

in LO and 

Di K ° +lf > NLO \x,Ql) = 0.50 a;- - 781 (1 



a;) 1 ' 58 , 



Df 0+K °> NLO \x,Ql)=Dl K0+K0 > NLO \x,Q 2 ) = 1.25a;- - 564 (1 - a;) 3 - 33 , 
Df° +r > NLO \x,Ql) = 0.99x-°- 601 (1 



aO 3 - 80 , 



D (K0+K ,NLO) (X; g2) = Q 33 ^-0.351 (1 _ ^0.65 (3) 



d (ko+k ,nlo) (X) g 2 } = 0.53 3,-0.571 (1 _ x y.m^ 

151 

in NLO. Here, it is understood that the values refer to the individual starting points 
given above. For the data that we fitted to, we find very small Xd.o.f. varues ; namely 0.53 
(0.52) at NLO (LO). The Xd.o.f. values achieved for the various data sets may be seen 
from Table I. Our FF also give a good description of the Z-resonance data from DELPHI 
and OPAL [|J, which yield just slightly larger values of Xd.o.f. The same is true for the 
lower-energy data taken by CELLO f2| and TASSO [|| at PETRA (^ = 35 GeV) and 



for the data collected by HRS p§ and TPC W§ at PEP (y/s = 29 GeV). Among the data 



that we compared with, those from CLEO [25] and ARGUS 0] have the lowest energy 



(y/s = 10 GeV). Only the ARGUS data give an exceptionally large Xd.o.f. 5 °f order 5. 

For the reader's convenience, we list simple parameterizations of the x and Q 2 depen- 
dences of our K° FF sets in the Appendix. We believe that such parameterizations are 
indispensable for practical purposes, especially at NLO. However, we should caution the 
reader that these parameterizations describe the evolution of the FF only approximately. 
Deviations in excess of 8% may occur for x < 0.1 and for Q > 100 GeV, in particular 
for the gluon. While this kind of accuracy is fully satisfactory for most applications, it 
is insufficient for the comparison with the high-statistics data collected at LEP. We wish 
to point out that all Xd.o.f. values presented in this paper have been computed using FF 
with explicit Q 2 evolution, which have an estimated relative error of less than 0.4%. 
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TABLE I. CM energies, experimental collaborations, numbers of data points used, and 
Xdof vames obtained at NLO and LO for the various e + e~ data samples discussed in the 
text. The data used in the fits are marked by an asterisk. 
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1.15 
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m 


4 


5.42 


5.32 



Since we have built in the cc and bb thresholds, we have three different starting scales 
Qo. To illustrate the relative size of the FF for the different quark flavours and the gluon, 
we have plotted them in Fig. |l] as functions of x for Q — 10 GeV. We show only the 
NLO results. The pattern is somewhat unusual and, contrary to naive expectations, not 
very similar to the K ± FF in our earlier work The u-quark, 6-quark, and gluon FF 
are rather hard, while the <i/s-quark and the c-quark distributions are soft. This pattern 
is already visible at the starting scales Q in Eqs. (fj) and (0), where we must keep in 
mind, however, that the starting scale Q takes on three distinct values for the light, c, 
and b quarks. Guided by our findings in connection with fragmentation, we would 
expect that, in Fig. p], the <i/s-quark FF should be hardest, and that the 6-quark FF 
should resemble that of the c quark. At this stage, it cannot be excluded that the relative 
importance of the individual quark flavours will need some adjustment. But for this we 
would need additional e + e~ data on inclusive K° production for which the fragmentation 
of the various quark flavours and the gluon is disentangled, similarly to what has been 
done in the case of charged-hadron production. Unfortunately, the existing information 
from ep and pp collisions does not help us much either. Due to its high threshold, b- 
quark production is absent at small pr, below 9.5 GeV. In our ep analysis, c/c production 
accounts for 8% (10%) of the cross section at pr = 5 (8) GeV, while, in our pp calculation 
for y^i = 1.8 TeV, its contribution at the same px values is 0.7% (0.8%), i.e., in both 
reactions it is small or negligible. 

The goodness of our fits to the ALEPH and MARK II |J data may be judged from 
Fig. |. At NLO (LO), we find xl.o.i. values of 0.60 (0.57) for ALEPH and 0.48 (0.49) for 
MARK II. 

The factorization theorem guarantees that the FF which we extracted from e + e~ 
data may also be used to predict other types of inclusive single K° production cross 
sections, e.g., for 77, ep, or hadron-hadron collisions. In the following, we shall present 
NLO predictions for inclusive photoproduction of Kg mesons at HERA and confront 
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them with preliminary data taken by HI [12|]. As in the HI measurement, we shall 
consider the spectrum of the produced Kg mesons, averaged over the rapidity range 
\y\ab\ < 1-5. We shall work at NLO in the MS scheme with Nf = 5 quark flavours, fix 
the renormalization and factorization scales by setting /i = M 7 = M p = = and 
adopt the NLO parton distribution functions (PDF) of the photon and the proton from 
Refs. |27|] and | 28| , respectively, together with our NLO FF. We wish to emphasize that 



also the hard-scattering cross sections will be calculated up to NLO. We shall evaluate a s 



to two loops with = 158 MeV [28|. The quasi-real photon spectrum will be simulated 



according to HI conditions, by imposing the cut 0.3 < z < 0.7 on z = E 1 / E e and choosing 
Qmax = 0-01 GeV 2 . Our predictions for £ = 1/2, 1, and 2 are confronted with the HI 
points in Fig. || The agreement is satisfactory as for both shape and normalization. 
Unfortunately, the HI data are accumulated at rather small pr {pr < 3 GeV), whereas 
our predictions should be more reliable at larger pt- We must keep in mind, however, 
that this represents the first measurement of inclusive K s production at HERA, based on 
data taken in 1993, and that the numbers are still preliminary. More data at larger p? is 
expected to appear after the analysis of the 1994 run is completed. As we see in Fig. |3|, 
the cross section shows only moderate scale dependence, which indicates relatively good 
perturbative stability. Notice that our prediction in Fig. |3] refers to K s production, which 
corresponds to the average of K° and K° . 

There only exists rather limited experimental information on inclusive Kg production 
in pp collisions. The only high-energy data available come from the UA5 Collaboration 



]TJ at the CERN SppS Collider and from the CDF Collaboration jTJ at the Fermilab 
Tevatron. In Fig. |], we show our predictions for the pt spectrum of pp — > Kg + X at 
y/s = 200, 546, and 900 GeV, with rapidity averaged over the interval —2.5 < y < 2.5. 
The calculation is performed at NLO in the MS scheme with Nf = 5 quark flavours 



using the CTEQ3 proton PDF [281]. The renormalization and fragmentation scales are 



identified and set equal to pr/2, pr, and 2pT- The agreement with the UA5 data [|13 
is quite satisfactory. It is best for the highest CM energy. Unfortunately, this data is 
accumulated at rather small px- On average, the agreement is best with scales equal to 
Pt- The data from CDF [14]] is more recent. It was taken at \fs = 630 and 1800 GeV. This 
data together with our theoretical results for scales pr/2, pr, and 2px are plotted versus 
Pt in Fig. |j| The experimental and theoretical results are both averaged over \y\ < 1.0. 
Unfortunately, the CDF data has rather small px, too. Again, the agreement of our 
calculation with the data is best for scales equal to pt- Only at very small pt, deviations 
from the data taken at \/s = 1800 GeV are noticeable. This is to be expected, since the 
theoretical predictions are valid only for large pt', their reliability at pt below 1.5 GeV 
is certainly questionable. At this point, we would like to encourage our experimental 
colleagues in the CDF Collaboration to also analyze the vast amount of data collected 
after 1989 with respect to light-meson fragmentation. In view of the considerable recent 
theoretical progress in this field, this would be interesting and exciting in its own right, 
rather than but a boring measure to assess backgrounds for certain other processes which 
presently happen to be more en vogue. In fact, this would allow us to test the QCD- 
improved parton model and, in particular, the factorization theorem at the quantum 
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level. 

Due to their limited px range and their modest accuracy, the data sets presented in 
Figs. |] and |5| are not so well suited for constraining the FF obtained from the e + e~ analysis. 
However, they provide a welcome cross check, in particular with respect to the gluon FF, 
which is only feebly constrained by the e + e~ data. To elaborate this point, we investigate 
the influence of the gluon fragmentation on the ep and pp cross sections. To that end, 
we repeat the calculations of Figs. |]-|5| switching off the quark FF. In Fig. |6], we show 
the outcome normalized to the full calculations for the ep cross section and the 200 GeV, 
630 GeV, and 1800 GeV pp cross sections. We observe that, in the \ow-px range, the pp 
cross sections are overwhelmingly dominated by the gluon FF. The ratio increases with 
CM energy and reaches 90% at the largest energy. This shows that, if it were not for 
the large errors, the pp data would be perfectly well suited for constraining the gluon 
FF. Looking back at Figs. [| and [5], it is fair to say that the strength of the gluon FF as 
obtained from our e + e _ — > /r* 1 + X fits is large enough to account for the pp data. This 
is in accord with recent studies of inclusive charged-hadron production in pp collisions 
2T| . We also examined in which x range the gluon FF maximally contributes to the pp 



inclusive cross sections in the considered pt range. Depending on the CM energy, the 
most important x values are concentrated around x = 0.4. This means that the pp data 
only constrains the gluon FF in a limited range of x. On the other hand, we know that 
the e + e~ data does not determine the gluon FF very accurately, i.e., a good description 
of the e + e~ data may also be obtained with a weaker gluon FF. The ep data also need 
a sufficiently strong gluon FF, in particular to describe the data near px = 2 GeV. At 
larger px, the influence of the gluon FF diminishes, and the quark FF come into play 
much more strongly. This is to be expected, since, in the ep cross section, the qg — > qg 
channel is similarly important as the gg — > gg and qg — > gq channels, even at small px- 

Having established the importance of the gluon FF for K s production in pp collisions, 
we should take a closer look at our assumptions concerning the gluon FF. These were 
twofold. We explicitly stated that we were going to assume (x, Qq) = (x, Qq) 
(a). A second, hidden assumption was that Df* had been well constrained by our previous 
analysis although only experimental information on the gluon FF for the sum of the 
charged hadrons had been available (b). By investigating the ratio of the cross section for 
inclusive kaon production in hadron collisions to that for charged hadrons, we may check 
both assumptions. For one thing, this ratio is approximately equal to the ratio of the 
respective gluon FF (at least for high CM energies), which enables us to test (b). On the 
other hand, this ratio has been measured for both charged and neutral kaons, providing 
us with a check of (a). In Fig. 0, we confront our predictions, based on assumptions (a) 
and (b), with the experimental data on neutral- kaon production in 1.8 TeV pp collisions 
by CDF [JyJ] and on charged-kaon production in 53 GeV and 27 GeV pp scattering by the 
British-Scandinavian Collaboration [30] at the CERN ISR and by the Chicago-Princeton 
Collaboration |31j at Fermilab, respectively. In the theoretical calculation of charged- 
hadron production, only charged pions and kaons are included. Protons, A hyperons, 
and other heavy hadrons are known to contribute little to the cross section and are 
neglected here. We find throughout good agreement with the data. All data, as well as 
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our predictions, approach a plateau at not-too-small pr- Its height is about 0.2, fairly 
independently of the CM energy or whether neutral or charged kaons are considered. 

At this point, we should compare our results on the K° FF with those obtained in 
Ref. [15|. The cross section of inclusive Kg production at y/s = 1.8 TeV under CDF 



conditions as predicted by Ref. [15] agrees reasonably well with our own calculation; it 



ranges between 50% and 70% of our result for p T between 2 and 7 GeV. 

These tests reassure us of the soundness of the assumptions concerning the gluon FF of 
the neutral kaons which we had to make in view of shortcomings in the presently available 
experimental information. From a theoretical point of view, it would certainly be desirable 
to constrain the K° FF by using just e + e~ data, as this would enable us to test them in 
other types of processes so as to probe the factorization theorem. Unfortunately, this is 
not yet possible, which has led us to use additional input to obtain FF that satisfactorily 
describe a variety of e + e~ , ep, and pp data. 



3 Summary and Conclusions 

We presented FF for neutral kaons, both at LO and NLO. They were constructed from 
fits to data on inclusive K° + K° production in e + e~ annihilation taken by MARK II |J 
at PEP (y/s = 29 GeV) and by ALEPH at LEP @. Although our FF were only fitted to 
the MARK II and ALEPH data, it turned out that they lead to an excellent description 
of other e + e~ data on inclusive K° + K production ranging from = 10 GeV to 
LEP energy. We always obtained Xdof values of order unity. The only exception, with 
Xd.o.f. ~ 5, occurred for the ARGUS data at y/s = 9.98 GeV. 

Since the e + e~ data do not constrain the gluon FF so well, we made NLO predictions 
for the pt spectra of Kg mesons produced inclusively in the scattering of quasi-real photons 
on protons under HERA conditions and in proton-antiproton collisions under UA5 and 
CDF conditions, and confronted them with the respective data. The agreement turned 
out to be quite satisfactory. We discovered that the gluon FF is very important to account 
for the pp data. We are thus faced with the unfortunate situation that the pp data almost 
exclusively tests the gluon FF, which is of little relevance for existing e + e~ data. Vice 
versa, the quark FF, which — up to a residual uncertainty in the relative importance of 
the individual flavours — are fixed by a wealth of e + e~ data, have hardly any impact on 
pp — > K s + X. The situation will be ameliorated as soon as the ep scattering experiments 
at HERA provide us with higher-statistics data, in particular at larger pr- In conclusion, 
present data does not yet allow us to test the universality of the FF in inclusive K° 
production; the situation rather requires that we exploit the universality postulated by 
the factorization theorem in order to extract meaningful FF. This was achieved in the 
work presented here. 

In order to make further progress, one would need e + e~ data on inclusive K° produc- 
tion in which the different quark flavours are tagged. Also, the gluon FF would have to be 
constrained better, e.g., by studying inclusive K° production in tagged three-jet events 
or by measuring the longitudinal part of the cross section, similarly to what has been 
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done for charged particles. As for ep and pp collisions, data at larger px with sufficient 
accuracy would be highly welcome, since they would allow us to quantitatively test the 
factorization theorem of fragmentation in the QCD-improved parton model. 
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A Parameterizations 

For the reader's convenience, we shall present here simple parameterizations of the x and 
Q 2 dependence of our FF.Q As usual, we introduce the scaling variable 

For A we use the MS value appropriate to Nf = 5 flavours, since the parameterization 

(5) 

would not benefit from the incorporation of discontinuities in s. is taken from our 
previous fit Q to be 107 MeV (195 MeV) in LO (NLO). Similarly to Eqs. we use 

three different values for Q , namely 

{V2 GeV, if a = u,d, s, g 

m(r] c ) = 2.9788 GeV, if a = c . (5) 

m(T) = 9.46037 GeV, if a = b 

This leads to three different definitions of s. For definiteness, we use the symbol s c for 
charm and Sb for bottom along with s for the residual partons. 

We parameterize our FF by simple functions in x with coefficients which we write as 
polynomials in s, s c , and Sb- We find that the template 

D(x,Q 2 ) = Nx a {l -xf (6) 



is sufficiently flexible, except for d( k °+ k : nl o) ? w here we include an additional factor (1 + 
7/2) on the right-hand side of Eq. (^). For s = s c = Sb = 0, the parameterizations agree 
with the respective ansatze in Eqs. The charm and bottom parameterizations 

must be put to zero by hand for s c < and Sb < 0, respectively. 

We list below the parameters to be inserted in Eq. (^) both at LO and NLO. The 
resulting parameterizations correctly describe the Q 2 evolution up to 8% for Q Q < Q < 
100 GeV and O.K x < 0.8. 



2 A FORTRAN subroutine that returns the FF for given x and Q 2 may be obtained from the authors 
via e-mail (binnewie@ipsl07.desy.de, kniehl@vms.mppmu.mpg.de). 
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LO FF for (K° + TC): 

• L>f 0+ ^°' L °)(x,g 2 ): 

N = 0.510 — 0.251s + 0.036s 2 
a = -0.841 - 0.285s + 0.021s 2 
= 1.550 + 0.712s - 0.069s 2 + 0.037s 3 

• ^f° + ^ LO) (x,g 2 ) = Df°+ s °."»(x,Q 2 ): 

iV= 1.470 - 1.088s + 0.276s 2 
a = -0.691 - 0.312s - 0.045s 2 + 0.038s 3 
/3 = 3.490 + 0.851s - 0.181s 2 + 0.098s 3 



D (K° + K°,LO)( xQ 2y 



N = 1.000 - 0.679s c + 0.111s 2 + 0.058s 3 
a = -0.738 - 0.302s c - 0.073s 2 + 0.084s 3 
P = 2.930 + 0.758s c - 0.118s 2 + 0.107s 3 



Dl K0+K °> u >\x,Cf): 



N = 0.680 - 0.470s b + 0.241s 2 - 0.115s 3 
a = -0.598 - 0.407s 6 + 0.059s 2 + 0.063s^ 
P = 1.930 + 0.554s 6 + 0.232s 2 - 0.186s 3 

• Df°+^> LO \x,Q 2 ): 

N = 0.430 - 0.881s + 0.860s 2 - 0.320s 3 
a = -0.374 - 2. 147s + 1.239s 2 - 0.232s 3 
P = 2.690 + 1.515s - 0.188s 2 - 0.071s 3 

NLO FF for (K° + K°): 

• L>f° + ^°' NL °)(x,g 2 ): 

N = 0.500 - 0.125s - 0.051s 2 
a = -0.781 - 0.500s + 0.264s 2 - 0.133s 3 
P = 1.580 + 1.074s - 0.380s 2 + 0.084s 3 
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Df ° +K °' NL0) (x,g 2 ) = Df +^ NLO )(*,Q 2 ): 



JV = 1.250 - 1.385s + 0.896s 2 - 0.251s 3 
a = -0.564 - 0.718s + 0.208s 2 

P = 3.330 + 0.100s + 0.388s 2 - 0.064s 3 (13) 



• Df° +K - NL °)(x,Q 2 ): 

N = 0.990 - 1.020s c + 0.639s 2 - 0.199s 3 
a = -0.601 - 0.666s c + 0.189s 2 

P = 3.800 + 0.192s c + 0.382s 2 - 0.077s 3 (14) 

• Dl K ° +lf > NLO \x,Qiy. 

N = 0.530 - 0.541s 6 + 0.429s 2 , - 0.171s 3 
a = -0.571 - 0.803s b + 0.371s 2 , 

(3 = 1.980 + 0.308s b + 0.298s^ - 0.035sj| (15) 

• Df 0+ *°> NL °)(x,Q 2 ): 

N = 0.330 - 0.307s + 0.075s 2 
a = -0.351 - 0.040s - 0.282s 2 + 0.033s 3 
= 0.650 + 2.141s - 0.457s 2 + 0.075s 3 
7 = 1.096s -0.198s 2 (16) 
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FIGURE CAPTIONS 



Figure 1: x dependence of the NLO set of K° + TC FF at Q 2 = 100 GeV 2 . 

Figure 2: Differential cross sections of inclusive K° + T( production at LO (dashed lines) 
and NLO (solid lines) as functions of x at ^Js = 91.2 and 29.0 GeV. The theoretical calcu- 
lations are compared with the respective experimental data by ALEPH B and MARK II 
||. For better separation, the distributions at 29.0 GeV have been divided by 10. 

Figure 3: The (preliminary) px spectrum of inclusive Kg production in ep collisions 
as measured by HI |12] is compared with the NLO calculation in the MS scheme with 



Nf = 5 flavours using the photon and proton PDF of Refs. |2yJ and |28[, respectively, 
together with our FF. The dashed/solid/dash-dotted curves correspond to the choices 
f = 0.5/1/2. 

Figure 4: The px spectra of inclusive Kg production in pp collisions as measured by 



UA5 |L3] at ^fs = 900, 546, and 200 GeV are compared with the respective NLO calcu- 
lations in the MS scheme with Nf = 5 flavours using the proton and antiproton PDF of 
Ref. P8|| . For better separation, the spectra have been separated by factors of 10. The 
dashed/solid/dash-dotted curves correspond to the choices £ = 0.5/1/2. 

Figure 5: Same as Fig. § but for data from CDF @] at y/s = 1800 and 630 GeV. 

Figure 6: Percentage of events with gluon fragmentation in the px spectra of K s mesons 
inclusively produced in pp and ep collisions. The solid, dash-dotted, dashed, and dotted 



lines represent our predictions for the 1800 GeV and 630 GeV CDF [jj, 200 GeV UA5 
l3| , and HI [12| experiments, respectively. 



Figure 7: Ratio of the differential cross section for inclusive kaon production to that for 
charged hadrons as a function of px- We compare the pp data on K§ mesons by CDF [fDH 
and the pp data on K + and K~ mesons (averaged) by the British- Scandinavian (BS) [[JIJ 
and Chicago-Princeton (CP) [|3T| Collaborations with the respective NLO calculations 
using our FF. We compute the denominators by summing over the charged pions and 
kaons. 
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